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Abstract 

The spin-1/2 zig-zag Heisenberg ladder (Ji — J2 model) is considered. A new representation 
for the model is found and a saddle point approximation over the spin-liquid order parameter 
(a„_i((Tn X ffn+i)) is performed. Corresponding effective action is derived and analytically 
analyzed. We observe the presence of phase transitions at values J2/ Ji = 0.231 and J2/J1 = 1/2. 
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1 Introduction 



Unconventional spin-liquid phases in frustrated spin chains attracted notable theoretical and 
experimental interest in recent years. The natural question is whether the frustrations in antifer- 
romagnetic Heisenberg chains can stabilize the new phases with exotic spin excitations observed 
in ladder systems. The model used to analyze the effect of frustration in antiferromagnetic spin 
chains is the so called spin-1/2 Ji — J2 model with the Hamiltonian 

H = Ji ^[(T„(T„+i - 1] + J2 y^[(7n'7ra+2 - 1], (1) 

n n 

where an = 2Sn are Pauli matrices. The bosonization analysis of this model was performed by 
Haldane [I], and the phase diagram has been studied intensively by various authors [21 [3l HJ EJ 
m [7] (for a review see also [8] ) . 

The interest in this model is not purely theoretical. There are inorganic compounds in nature, 
such as CS2CUCI4 L9J, CuGeOs [10], LiVaOs [H], or SrCu02 [12], which can be described by 
the spin-1/2 Ji — J2 chain Hamiltonian Eq. ([T]). 

The investigation of the spin-1/2 Ji — J2 model's phase diagram ([1]) starts from the weak 
coupling limit when J2 « Ji- Classically, in this limit, the Ji — J2 chain has a long range 
Neel order. The excitations are massless spin-waves, frustrated by an irrelevant perturbation. 
At the values J2/J1 > 1/4, the spins in the ground state are arranged in a canted configuration, 
in which each spin makes a fixed angle a = arccos[— J1/4J2] with its predecessor. The classical 
ground state of the model is doubly degenerate since the spin configurations can turn clockwise 
and counterclockwise with the same energy. 

In quantum field theory, it is believed that at larger values, J2/J1 > 1/4, a phase transition 
of the Berezinskii-Kosterlitz-Thouless (BKT) type [13] occurs, which separates the gapless spin- 
1/2 Heisenberg phase from a fully massive region. This phase is characterized by the two- 
fold degenerate dimerized ground state and a spontaneous breaking of the lattice translation 
symmetry. Frustration stabilizes this gapful phase. The actual value for the ratio of the coupling 
constants at the transition point was found numerically in [21 [Sj HI EJ [6] to be slightly lower than 
J2/J1 = 0.241 due to quantum fluctuations. 

At the larger value (when J2 = 1/2 Ji > 0) the model coincides with the Majumdar-Ghosh 
(MG) model [HI [l5] and is exactly integrable [Ml (ITj. The existence of a mass gap in the 
MG model has been shown rigorously in [18]. The correlation function is found to be zero at 
distances larger than the lattice spacing. The ground state of the phase at J2/J1 > 1/2 is found 
to be a condensate of dimerized singlets of pairs of neighbor spins, which is Z2 degenerate. This 
spontaneous breaking of the Z2 discrete symmetry creates a kink type of topological excitation, 
the tails of which end in the different Z2 vacua at ±00. The spectrum of this excitation is massive, 
and the gap is decreasing at J2/J1 00, when the system becomes a pair of noninteracting 
spin - 1/2 Heisenberg chains. 

In the present article, we develop an approach to analyze the model on the matter of critical 
behavior, based on the idea that the middle phase, 1/4 < J2/J1 < 1/2, can be characterized by 
the spin-liquid order parameter 

V = {^n-li^n X ^n+l)) (2) 

defined on the triangles of the zig-zag chain. Performing a mean-field (saddle point) approx- 
imation, we reduce the model to an extended Heisenberg chain with topological term, which 
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appears to be integrable. Then, by use of the technique of Thermodynamic Bethe Ansatz, we 
calculate the effective action (formula (j34p ) as a function of the spin-liquid order parameter ^p. 
The analysis of the effective action shows that, at the points J2/J1 = 0.230971 and J2/ J\ = 1/2, 
we indeed have phase transitions, as it was expected. 



2 New representation for the Ji — J2 model 

For an alternative treatment of the Hamiltonian ([1]), we make use of the following identity 
for the square of the Hermitian operator Xahc = ^ai^b x <?c) (scalar chirality operator) 



{aa{ab X ac)f = -2[[aaSb - 1] + [^a^c - 1] + [^b^c - 1] j ■ (3) 

On the right hand side of this identity, one may recognize a sum of three different spin exchange 
terms, marked by mutually non-coinciding indices a ^ h ^ c. The form of Eq. ^ suggests 
a transformation for the Ji — J2 Hamiltonian, a dual representation of which would contain 
the square of the above mentioned scalar chirality operator. Thus, using Eq. ([3|), we map the 
expression ([1]) to 

H = Yy^n^n+l - 1] - I X] ^rln+l,n+2^ 9 = j "^^j , (4) 

n n 12 

were we have rescaled the Hamiltonian H by the constant factor (Ji — 2J2), i.e.. 

This will be our starting point. We would like to emphasize however, that, to the best of our 
knowledge, the Heisenberg zigzag ladder ( Ji — J2 model) has not been represented in this form 
in the literature previously. 



3 Mean-field theory and its Bethe Ansatz solution 

Our further analysis is close in spirit to that of Affleck and Marston, applied in Ref. [19] for 
the solution of the Heisenberg-Hubbard model. It is based on the approximation of the partition 
function of the given model by an exactly integrable one. Namely, we introduce a scalar field, 
(pn, by a Hubbard-Stratanovich transformation, which maps the Hamiltonian H to Ti., where 



= H + - "^{(pn - Xn,n+l,n+2f = ^[(Tn-Tn+l " 1] " fi' X] ^nXn,n+l,n+2 + 2 X] '^"^ 
n n n n 

This map induces a constant factor in the partition function 

Z = Tr (exp = const Tr (^j exp ({(/>„}) JJ d</.„^ , (7) 

which however is irrelevant. The fact that the transformation ^ leaves the dynamics of the 
model unchanged may be shown using the coherent state path integral representation for the 
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partition function, where the functional integral over the field {cpn} can be exactly evaluated. 
Below, we will investigate the zero temperature limit, /3 — > oo, of the integral on the right hand 
side of Eq. ([7|) in the saddle point approximation. We are going to analyze the mean-field theory 
corresponding to a certain saddle point, which we believe gives the main contribution. The 
question of the existence of any other saddle points, however, will be left for future investigations. 
More precisely, we consider the saddle point equation 

The solution of Eq. ^ with regard to the bosonic fields {(pn} can be obtained by the substitution 
of Eq. ([6]) into Eq. dH) . It has the form of a set of (the lattice size) coupled equations 

Tr(Y U.1 _j.oe-'3^({'^"»'l 

- Tr (e-/3W({*'^})) " (Xn,n+l,n+2), (9) 

where n = 1 . . . A^, and we have cyclic boundary conditions. It would be reasonable to restrict 
ourselves by some (quasi) translational invariant, homogeneous saddle points. Therefore, we 
consider the solution where the operators Xn,n+i,n+2i for all n = 1 . . . A^, have the same mean 
value. Then the set of coupled equations ([9]) simplifies, and for all values of n acquires the form 

{Xn,n+l,n+2) = (10) 

Thus, in this way, the original problem reduces to the eigenfunction problem for the mean-field 
Hamiltonian 

= ^[o'nO'n+l - 1] - a^'?n(<?n+l X o'n+2), a = Qlf. (11) 
n n 

This model appears to be exactly solvable by means of Bethe Ansatz. This is because the 
second term in expression (jlip commutes with the first term, which, in turn, is the Heisenberg 
Hamiltonian. Therefore, skipping the demonstration of the standard technique of Algebraic 
Bethe Ansatz (since it repeats the one for the XX Z model |17j). we present here only the 
solution. The eigenvectors can be parameterized through the set of parameters (rapidities), 
{xi}, which satisfy the set of Bethe equations 

_ -xN M _ _ 

Xj + i J Xj - Xk + 2i' 

The corresponding state has a total spin projection Sz = N — M and energy 

M g 

E{xu ..,xm) = (1 + 2a9,,) r- (13) 
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4 Free energy 



In the present section, we present our calculations of the ground state energy of the Ji — J2 
model in the mean-field approximation. Namely, we calculate the free energy in the model 
The model (jlip has been studied in Ref. [20] for a fixed value of the parameter a. Here, we 
present our exact analytical calculations of the free energy of the mean- field model and 
analyze, in detail, the free energy as a function of the parameter a. The calculation is based on 
the method of Thermodynamic Bethe Ansatz, introduced in (for details see also [22J). By 
definition, the thermodynamic limit is given by the following conditions 

M 

N ^ 00, M ^ 00, ^ = const. (14) 



In the thermodynamic limit (jl4p . the Bethe equations ()12p become integral equations. In order 
to represent these integral equations in a convenient form, we introduce the densities p{t), defined 

as p{t) = (^^\t=t{x)j ■ Then the equations (fT2]) reduce to the following system of N integral 
equations for the densities 

Unix) = Pn{x) + Pn{x) + ^^Tjk * Pk{x). (15) 

k 

Here we introduced the notations 

In 

an[x) = — 2, (16) 

for n = I...N. The functions pn{x) and Pnix), which are unknown, denote particle and 
hole densities, respectively. The index n represents their correspondence to n-strings. The 
convolution operation, " * " , is defined as 

/ + 00 
fix - y)g{y)dy, (17) 
'OO 

for any given pair of functions / and g. This is the conventional definition. The functions Tnm{x) 
for n, m = 1 . . . N , involve the expressions an{x) in their definition. They have the following 
form 



-^nm(-^) — * 



' a\n-m\{x) + 2a\n-m\+2{x) + 2a\n_ra\+i{x) + ... 

... + 2an+m-2{x) + an+m{x) for TL ^ 771, ^^^^ 

2a2{x) + 2a4{x) + ... + 2a2n-2{x) + a2n{x) for n = m. 



As we have already mentioned, the equations (jlSp represent the thermodynamic limit of 
the Bethe Equations. These are integral equations with respect to particle and hole densities, 
containing all of the information about the energy spectrum. Suppose that the system is in a 
state characterized by densities Pj{x) and Pj{x). Then the equilibrium dynamics of the system 
at temperature T can be extracted by minimizing the free energy, F = E — TS, with respect 
to the independent pj. This yields the following non-linear integral equations for functions 

r/„(x) = pn{x)/pn{x), 
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oo 

lnr/„ = ^ + J]]r„fc*ln(l + 7?^^), 5„ = -Svr (1 + 2a9^) a^. (19) 
fc=i 

In order to analyze the ground state energy of the mean-field model Eq. (jlip which is under 
our current consideration, we need to go to the zero temperature limit in Eq. (|19p . For this 
purpose, let us introduce a set of new functions, en(x), as r]n{x) = exp{e„(a;)/r}, n = 1 . . . A^, 
and substitute them into Eq. (I19p . Then, in the zero temperature limit, T — > 0, Eq. (I19p acquires 
the following form 

e^{x) = -Svr (1 + 2adx) s{x) + s* e^(x), (20) 
e„(2;) = s * (e„„i + e„+i)(x), n > 2, 

where the function s(x) is defined as s(x) = [4 cosh(7rx/2)]~^. The action of the dagger (minus), 
t (-), in Eq. ([20p . leaves only the positive (negative) part of the corresponding function, e„(x), 
as 

en{x) > 
e„(x) < 0, 

(21) 

By definition, all e„(x) are measured in units of kT (where we set A; = 1), and therefore have 
magnitudes of energy. Equations (l2Qp unambiguously define the solutions for functions 6^ (x) 
provided that 



e„(x) if 

if 

e~{x) = e„(x) - ei{x). 



lim = 2B, (22) 

n—*oo n 

where B is the "magnetic field", which in our case, Eq. (jlip . is zero. It is transparent from 
Eq. (j20|) . that e„(x) > for n = 2,3,.., and only the function ei(x) can be positive, as well as 
negative (can change its sign crossing the x axis at a certain point). The solution of the system 
(j20p can be then expressed in terms of ei(x) as 



e„(x) = ei{x) = a„_i * e{{x) + 2(n - 1)5, n = 2, 3... 



ei(x) = -Svr (1 + 2a9^) s(x) + / {s * ai){x - y)ei{y)dy. 

Jei>0 

From the last equation we see that if there exists such a point, x = a, where the function ei(x) 
changes its sign (and, therefore, ei(a) = 0), then ira > 1. Thus, one arrives at a Wiener-Hopf 
type integral equation for the unknown function ei(x) 

ei(x)=eo(x)+ j R{x-y)ei{y)dy, (23) 
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where 



r . p-M 

eo{x) = -87T{l + 2ad,)s{x), R{x) = / — e"^"^- — . 

J 27: 2 cosh uj 

—oo 

The same kind of non-linear integral equation, occurring when the so called disturbance term, 
eo(x), is not an even function of x and changes its sign, appears in the solutions of staggered 
zigzag ladders with broken one-step translation symmetry ^23j, [24j. Therefore, according to our 
experience drawn from the previous works, we assume that €i{x) = for x < a. Then, Eq. (|23p 
will be valid for x > a. This assumption does not affect the right hand side of Eq. (j23p and, 
without loss of generality, gives the same solution. This solution can be obtained following the 
standard steps of the technique of Wiener-Hopf integral equations. First, we apply a Fourier 
transform 

oo oo 

fix) = I /» = I dxe--/(x), 

— oo — oo 

to the functions eo(a;), ei(x) and R{x). The substitution of these functions in the form of Fourier 
integrals into Eq. (p3]) yields 

oo 

/ ^e~*""{[l-i?M]eiM-eoM} =0, x>a. (24) 

— oo 

The equation (|24p can be rewritten in an equivalent form, as 

[1 - R{u:)\h{^) - eo(w) = e*""/i-(^), (25) 

where h±{uj) are the boundary values of analytic functions which do not have poles in the 
upper (n_|_) and lower (n_) complex half-planes respectively, and have zero limiting values at 
corresponding infinite points. Hence, with our assumption, we will have ei{ui) = e*'^"e+(w). The 
kernel in Eq. ()24p can be factorized. It is precisely this factorization property of the kernel which 
is responsible for the solvability of Eq. (|24p . This means that the kernel can be represented as 
a product 

[G,MG-Hl-'.i-iJM = ^^, (26) 

where G±{uj) are the boundary values of the analytic functions which do not have zeroes or 
poles on II-i-, respectively, and have the property G^{oo) = G^{oo) = 1. Then, for the Fourier 
components ei{io) in Eq. (j23p . the solution, when x > a, will be 

6+(u;) = G+(u;)P+ [G_(a;)e-*-'^eo(^)] • (27) 

Here, the operators P± are projectors, defined as P±[f{uj)] = f±{u!), for any given function 
f{x). For example, the action of the projector P+ on the sum of the Fourier components of the 
function f{x) and a constant c, yields 

/ + 00 "I Z' + OO 

dxe^'^^/(x) =c+ / dxe'^'^fix). (28) 
-oo J Jo 
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-daf/a 




0.2 0.4 0.6 0.8 

Figure 1: The function —daF{a) /a versus a 



The complex functions G+ (uj) and (uj), from the factorization equation (I26p , can be calculated 
exactly. They have the following algebraic forms 

^ ^ ^ e^p{-i^ +i^ln(i^)} , , , , 

Now one can derive the explicit form of the solution (|27p for a > 0. In order to do this, 
one just has to make use of the above mentioned property of projectors P+ and P_, given by 
Eq. (|28|) . Namely, upon application of Eq. (|28|) to the expression in brackets in the right hand 
side of Eq. (j27p . one will express the solution Eq. (j27p in the form of an infinite sum 



e+M = -ii.G^iu;) E(-l)'[l - + 1/2)]- ^^t{k + l/2) 



The first two terms of the sum in this equation have been obtained in Ref. [20]. Now, in order 
to calculate the free energy in the mean- field model Eq. (jlip . we need to find the parameter 
a = a{a), defined by the condition ei(a) = 0. This condition can be rewritten as 

/duj 
7r-e+(^) - ^ , lim (^e+{uj). (31) 

Substituting the solution Eq. (|3Up for e+(f-i^) into Eq. (j3ip . we will represent Eq. (|3ip in an 
equivalent form as 

oo 

J]](-l)^[l - 2TTa{k + l/2)]e^™('=+^/2)^_(_^^(^ ^ 1/2)) = 0. (32) 

k=0 

The free energy per site at T = 0, which is the ground state energy, will explicitly depend 
on a. From the definition of the free energy, we have 

F{a) = -(41og2-l)-/ s{x)ei{x)dx 

JxeL+ 

oo 

= -(41og2-l)- I ^e--'^6+(^).+ (-c.). (33) 
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9=0.6 

Figure 2: Vacuum energy Eq versus parameter ip for different values of g. 



Here, the integration range of the first integral, L+, is given by those values of x, where ei(x) is 
positive. Thus, substituting Eq. ([30|) into Eq. ([3T]) . we get the exact ground state energy of the 
model Eq. (fTTIl as 

oo 

F{a) = 4- 41og2 + 2 ^ - 2to(A: + 1/2)] (34) 

n,fc=0 

e-^"(^+"+i)G+(i7r(n + 1/2))G„ (-i7r(A; + 1/2)) 
k + n + l 

5 Effective action and phase transitions 

The expression for the energy per site, can be considered as an effective action for the 
order parameter 99, defined by Eq. ([2]). The explicit form of the effective action, Eq. (j34p . allows 
for further investigations; in particular, with regard to the matter of critical behavior, one can 
analyze in details the saddle point equation corresponding to our mean-field theory. In terms of 
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the parameter a, the saddle point equation ([9]) reads 



daF{a) 



a/g, 



or, if a 7^ 



daF{a) 



i/g. 



(35) 



a 



For any given g, this equation always has a zero solution a = 0. In order to find a non-zero 
solution, we performed a numerical evaluation of the function —daF{a)/a (where the analytical 
form of F{a) is given by Eq (j34p ). with the condition given by Eq. (|32p . The plot is presented 
in FigHJ From this picture, one can conclude that Eq. (j35|) has a solution when and only when 
g exceeds the critical value gc, where 1/gc equals to the maximal value of —daf{a)/a. For that 
value, our calculations give gc = 0.428646. Then, from Eq. (jj]), one can find the corresponding 
critical ratio 



which is in a good agreement with the expected value. It is also interesting to investigate the 
behavior of the effective potential, Eq. (|34|) versus the order parameter if for different values of 
g. The plots are presented in Figl2j When g is less than gc, there is only one vacuum energy 
minimum at 93 = while, for g > g^., two new minima appear. There exists another value of g, 
which we mark as gc2 {gc2 = 0.555083, [J2/Ji]c2 — 0.263052), which occur when the magnitudes 
of three vacuum energy minima are the same. At this point we do not have an additional phase 
transition, since the order parameter, ip, is smooth and finite. However, it would be interesting 
to understand the reason and consequences of such behavior. 

Two non-zero minima become infinitely deeper and the positions of minima approach to zero 
from both, left and right hand sides, as g further goes up to -|-oo. This scenario corresponds to 
the Majumdar- Ghosh limit, J2/J1 1/2, suggesting the next phase transition where the chiral 
order parameter ip vanishes and another, fully dimerised phase appears. Thus, our description of 
the intermediate phase with non-zero chiral order parameter ip complements the understanding 
of the fully dimerised phase for J2/J1 > 1/2. 

Concluding remarks. We have derived an effective action for the spin-1/2 Ji — J2 model as 
a function of the spin-liquid order parameter ip = (CTn-i(<?n x (^n+i)), and have observed the 
presence of two phase transitions at points close to the expected values: (i) when {J2/J1) < 
0.230791, we have an ordinary critical phase of isotropic Heisenberg model; (ii) in the middle 
phase, when 0.230791 < (J2/J1) < 1/2, the ground state is Z2 degenerate with two signs of 
order parameter. Due to this degeneracy, kink-like topological excitations are present and their 
condensation may characterize the third phase at J2/J1 > 1/2. Though the described picture 
do not coincide, but at the same time is not in contradiction with the well known description 
of this phase in the thermodynamic limit, when the two states, one with wave vector k = 
(ground state for finite system) and another with k = it (first exited state for finite system), 
collapse to each other and give rise to the dimerization pattern (two-fold degeneracy) and the 
breaking of one-step translational invariance. 

In our opinion, the developed approach based on the chiral order parameter ip, is alterna- 
tive to the known methods for description of the intermediate phase and provides promising 
possibility to investigate this important problem further. 
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